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ABSTRACT 



Purpose: The objective of this study was to assess the trabecular bone changes that occurred around functional 
endosseous dental implants by means of radiographic image analysis techniques. 

Materials and Methods: Immediate preoperative and postoperative periapical radiographs of de-identified implant 
patients at the University Dental Hospital of Manchester were retrieved, screened for specific inclusion criteria, 
digitized, and quantified for structural elements of the trabecular bone around the endosseous implants, by using 
image analysis techniques. Data were analyzed using SPSS version 11.5. P values of less than 0.05 were considered 
statistically significant. 

Results: A total of 12 implants from 1 1 patients were selected for the study, and 26 regions of interest were obtained. 
There was a significant increase in the bone area in terms of the mean distance between nodes (p =0.006) and a 
significant decrease in the marrow area in terms of the bone area (p =0.006) and the length of marrow spaces (p= 
0.032). 

Conclusion: It appeared that the bone around the implant underwent remodeling that resulted in a net increase in 
bone after implant placement. (Imaging Sci Dent 2014; 44: 129-35) 
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Introduction 

The success of dental implants can be assessed clinically 
and radiologically. In terms of the radiological assessment, 
the features of successful implantation include a mean 
crestal bone loss of no greater than 1 .6 mm during the first 
postoperative year and no greater than 0.13 mm through 
the maintenance period, and lack of peri-implant radiolu- 
cency. 1 ' 2 However, we have little knowledge of the effects 
produced by functional endosseous implants on the quality 
and architecture of the surrounding trabecular bone. 

Imaging is an essential part of endosseous implant ther- 
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apy in preoperative, intraoperative, and follow-up assess- 
ments. An ideal imaging technique should reveal many 
significant characteristics at specific stages of the implant 
therapy. Therefore, a combination of radiograph techni- 
ques is used with certain indications. 3 

Digital radiography is known to be one of the most sen- 
sitive methods for assessing subtle changes in trabecular 
bone density. It has high image photometric accuracy, wide 
dynamic range, and high resolution capacity for dental 
assignments. 3 Furthermore, it allows for image enhance- 
ment and processing on a computer monitor, such as mag- 
nification, skeletonization, and contrast optimization. 3 " 5 

The objective of this study was to assess the trabecular 
bone changes that occurred around functional endosseous 
dental implants by means of radiographic image analysis 
techniques. A specific image analysis technique was per- 
formed in order to quantify the structural elements of the 
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trabecular bone. Accordingly, the null hypothesis was as 
follows: There would be no significant quantifiable differ- 
ences in the trabecular bone structure on radiographs sub- 
jected to repeated analyses. 

Materials and Methods 

Selection of patients 

The computer database of anonymized implant patients at 
the University Dental Hospital of Manchester was retriev- 
ed and screened for the following inclusion criteria: 1) 
endosseous implants of the Frialit-2 implant system (Fri- 
adent GmbH, Mannheim, Germany), which was the most 
commonly used implant type in the institution, 2) avail- 
ability of (immediate) preoperative and (12-24 months) 
postoperative periapical radiographs, 3) adequate techni- 
cal quality of radiographs, 4) clinical and radiographic 
evidence of the absence of local pathology, and 5) patients 
with unremarkable past medical history. 

Image digitization 

The periapical radiographs were digitized by an Epson 
Expression 1680 Pro Scanner (Epson UK Ltd., Hemel 
Hempstead, UK) using an 8-bit gray scale and with a res- 
olution of 1200 dots per inch (dpi) so that they could be sub- 
jected to several software processing techniques, through 
a number of stages. In the first stage, Adobe Photoshop 
Elements version 3 (Adobe Systems, San Jose, CA, USA) 
was used to select three boxed regions of interest (ROIs) 
in each periapical image, first by using the marquee tool 
and then by cropping. For meeting the inclusion criteria of 
the ROI, one region was of 150 x 150 pixels in the tagged 
image file (TIF) format just below the root apex for single- 
rooted teeth and at a mid-point below and between the 
apices for multi-rooted teeth. The other two regions of 
150 (height) x75 (width) pixels in the TIF format were 
selected parallel to the apical third of the root on both 
mesial and distal aspects (Fig. 1). In general, the three 
boxes were confined to the periapical areas without lamina 
dura containing the trabecular bone seen as white lace-like 
structures on the gray background of the marrow. As the 
final step in this stage, the three cropped images were con- 
verted into an Adobe Photoshop grayscale image and then 
filtered (radius 31, amount 50, and threshold 0) (Fig. 2). 
Image filtration could be defined as a procedure in which 
a computer adjusts an image by convolving it with a filter 
function, in terms of smoothening and sharpening in order 
to improve the image quality. 3 




Fig. 1. Digitization, selection, and cropping of the images, (a) A 
periapical radiograph is digitized by a scanner using an 8-bit gray 
scale and a resolution of 1200 dpi. (b) A magnified region of inter- 
est (ROI) of 150 x 150 pixels in the tagged image file (TIF) format 
below the root apex is cropped by Adobe Photoshop Elements 
version 3. (c) A magnified ROI of 150 x 75 pixels in the TIF at the 
mesial aspect of the root apex, (d) A magnified ROI of 150 x75 
pixels in the TIF format at the distal aspect of the root apex. 

Image optimization 

In the second stage, the 30 filtered images were processed 
using EMAGO Software (Oral Diagnostic Systems, Ams- 
terdam, the Netherlands) in order to obtain contrast stretch- 
ing followed by optimization of the contrast of gray levels 
from which the histograms for each image were acquired. 
The peak of each histogram was identified manually to 
determine the threshold numbers for subsequent binariza- 
tion (i.e., segmenting the image into trabeculae and mar- 
row). 5 On some occasions, there were two peaks for one 
histogram, and in these situations, a subjective estimate 
was made to identify the midpoint between the peaks for 
the threshold number determination. Then, the images 
were binarized, and subsequently, the noise was suppressed 
(Fig. 2). 

Trabecular analysis system 

In the last stage, the Win TAS (S. Paxton, University of 
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Fig. 2. The images are the procedures for image digitization, optimization, and analysis (marrow). 



Leeds, UK) automated trabecular analysis system was used. 
The first step in this stage was to apply thresholding (i.e., 
separating the trabecular network from the background 
marrow) in order to further suppress the image noise so 
that they are recognized by the software. 6 As a result, at 
this level, the marrow was seen as white and the trabecu- 
lae as black. Then, the images were eroded (i.e., pixels 
from the edges of objects in binary images were removed 
by separating objects that touched each other and remov- 
ing isolated pixels) once and dilated (i.e., pixels at the 
edges of the objects were added by joining discontinuous 
objects and filling in the holes) once. 3 Finally, the result- 
ing images were analyzed and skeletonized (i.e., repeated 
elimination of pixels from the borders of objects in a bina- 
ry image until they are diminished to single pixel-wide 
skeletons) to their final form of algorithmic "parallel plate" 
parameters (Table 1). The gray scale of the images was 
inverted to make the trabecular component white, and thus, 
its measurements were obtained. The last stage was repeat- 



Table 1 . Parallel plate parameters or histomorphometric variables 
as listed in the data printout from Win TAS 



Abbreviations 



Name/definition 



BV 


Bone volume 


BS 


Bone surface 


TbTh 


Trabecular thickness 


NNd 


Node number: Number of branch point 




in a trabecular network that has been 




skeletonized. 


NTm 


Terminus number: Number of end point 




in a trabecular network that has been 




skeletonized. 


NNd:NTm 


Node terminus ratio: Index of spatial 




connectivity. 


BAr 


Bone area 


BPm 


Bone perimeter or boundary 


TbWi 


Trabecular width 


TbN 


Trabecular number 


TbSp 


Trabecular separation 


Tar 


Tissue area 


BS3d 


Bone surface in three dimensions 


SDn 


Surface density 


TbLe 


Trabecular length 
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Fig. 3. Inverted image is analyzed (trabecular) using Win TAS for the images shown in Fig. 2. 



Table 2. Cases selected for the study 



Number of study 


Age at the time 
of the first radiograph 
(years) 


Sex 


Site of implant 


Time between the first 

and the second 
radiographs (months) 


1 


23 


F 


Rt. Max. C and PM2 


21 


2 


43 


F 


Rt. Max. LI 


15 


3 


50 


M 


Rt. Max. CI 


20 


4 


34 


M 


Lt. Max. CI 


18 


5 


20 


M 


Rt. Max.CI 


19 


6 


62 


M 


Rt. Max. Ml 


16 


7 


20 


F 


Lt. Max. CI 


17 


8 


23 


M 


Rt. Max. C 


19 


9 


36 


F 


Rt. Max. C 


15 


10 


20 


M 


Rt. Max. PM2 


18 


11 


23 


M 


Lt. Max. PM1 


24 



F: female, M: male, Rt: right side, Lt: left side, Max: maxillary, CI: central incisor, LI: lateral incisor, C: canine, PM1: first premolar, PM2: second pre- 
molar, M 1 : first molar 



ed for the whole set of images (Fig. 3). As a result, two 
sets of skeletonized images led to two sets of quantitative 
tables (marrow and trabecular). 

Statistical analysis 

Data were analyzed in SPSS version 11.5 (SPSS Inc., 
Chicago, IL, USA). The mean value, standard deviation 
range, and 95% confidence interval were calculated for 
each parameter obtained by Win TAS (marrow and trabe- 
cular) for the first and the second sets of images. The pair- 
ed data were then analyzed to calculate the repeatability of 



parameters. Two-tailed paired t-tests were used to assess 
the significance of differences between the first and the 
second image analyses. The p values of less than 0.05 were 
accepted as statistically significant. 

Results 

A total of 12 implants from 11 patients were selected 
for the study; however, only 26 ROIs (instead of 36) were 
finally included because the space available for ROI place- 
ment was insufficient in the rest of the cases (Table 2). 
There was a significant increase in bone area in terms of 
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Table 3. A comparison between the mean "bone" data from the 
first (baseline) and the second (follow-up) radiographs (n=ll) 



Measurement 


Mean at baseline 


Mean at follow-up 


NNd 


29.4 


33.6 


NTm 


23.5 


24.6 


BAr 


5815 


7913.9 


BPm 


1076.9 


1206.9 


TbWi 


10651 


12903 


TbN 


0.047 


0.055 


TbSp 


19201 


9176 


BS3d 


0.066 


0.076 


SDn 


0.203 


0.147 


Mean Nd-Nd 


22.3 


31.6 


Mean Tm-Tm 


14.3 


13.7 


Mean Tm-Nd 


16.6 


20.2 


TbLe 


622.5 


753.5 



Mean at baseline, the mean of the first set of radiographs at the time of 
inserting the implant; mean at follow-up, the mean of the second set of 
radiographs at the time of making the implant functional; NNd, node 
number; NTm, terminus number; BAr, bone area; BPm, bone perimeter; 
TbWi, trabecular width; TbN, trabecular number; TbSp, trabecular separa- 
tion; BS3d, bone surface; SDn, surface density; mean Nd-Nd, mean dis- 
tance between nodes; mean Tm-Tm, mean distance between terminals; 
mean Tm-Nd, mean distance between terminal and node; TbLe, trabecular 
length 



Table 4. A comparison between the mean "marrow" data from the 
first (baseline) and the second (follow-up) radiographs (n=ll) 



Measurement 


Mean at baseline 


Mean at follow-up 


NNd 


41 


28.7 


NTm 


21.5 


23.9 


BAr 


6989 


5367 


BPm 


1169 


1091 


TbWi 


12687 


10219 


TbN 


0.054 


0.050 


TbSp 


10271 


15950 


BS3d 


0.075 


0.070 


SDn 


0.172 


0.211 


Mean Nd-Nd 


25.4 


21.2 


Mean Tm-Tm 


13.2 


13.5 


Mean Tm-Nd 


15.6 


16.1 


TbLe 


772 


573 



Mean at baseline: mean of the first set of radiographs at the time of insert- 
ing the implant. Mean at follow-up: mean of the second set of radiographs 
at the time of making the implant functional. NNd: node number. NTm: 
terminus number. BAr: bone area. BPm: bone perimeter. TbWi: trabecular 
width. TbN: trabecular number. TbSp: trabecular separation. BS3d: bone 
surface. SDn: surface density. Mean Nd-Nd: mean distance between nodes. 
Mean Tm-Tm: mean distance between terminals. Mean Tm-Nd: mean 
distance between terminal and node. TbLe: trabecular length 

the mean distance between nodes (p =0.006) and a signifi- 
cant decrease in marrow area in terms of the bone area 
(p=0.006) and the length of marrow spaces (p=0.032). 
Tables 3 and 4 show the mean bone and marrow parame- 
ters from the baseline (first) and follow-up (second) sets 
of radiographs for ease of comparison. 



Discussion 

The development of a non-invasive technique for mea- 
suring bone architectural changes has been a challenge 
for many research groups. 7 Texture analysis of peripheral 
microtomographs has been suggested as an alternative 
approach. 7 However, in a dental context, these techniques 
are difficult for many clinical research groups because of 
the limited availability of equipment or practical problems 
in their use. However, although two-dimensional radio- 
graphic measurements fail to spot some information regard- 
ing the out-of-plane curvature of trabeculae, they have 
significant correlation with data provided by three-dimen- 
sional computed microtomography. Furthermore, it has 
been shown that trabecular bone microarchitecture can be 
measured by the texture analysis of radiographs of the 
dried bone. 8 

There were some technical limitations in this clinical 
study that should be addressed. The first was that the 
images had to be modified in Adobe Photoshop, processed 
in EM AGO and then analyzed in Win TAS. This long 
course of action was rather time consuming and tedious 
in terms of filing and saving the images, and confusing in 
terms of transferring them from one system to another. 
This difficulty could have been overcome if a special soft- 
ware system had been designed specifically for dental 
radiology practice. 

The second limitation was the size and inflexibility of 
the rectangular marquee tool used in the Adobe Photoshop 
elements, which was larger than the implant and as a result, 
could not follow the implant contour. This led, in some 
cases, to the elimination of one or two out of the three 
ROIs (a, b, and c). However, the elimination was consis- 
tent through the first and the second sets of the image anal- 
ysis. A solution for this problem might be achieved by 
designing the marquee tool to be freely variable in size 
and shape, perhaps by using a cursor-generated freehand 
ROI. 

Third, although images were acquired by a paralleling 
technique with the same X-ray machine and fixed expo- 
sure parameters, there was still some variation in the X- 
ray angulation that resulted in a minor difficulty of retain- 
ing the same proportions of ROI on the paired images. In 
some cases, there was a shortage of space around the 
implant on one radiograph, such that it was difficult to 
obtain the same ROI, leading to the exclusion of that ROI 
from the analysis or a change in the position of the mar- 
quee tool in others. The importance of standardization was 
highlighted by Jonasson et al, particularly, as mentioned 
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above, in relation to the angulation of the X-ray beam. 9 
Grondahl et al have mentioned that even minor deviations 
in the vertical angulations of the X-ray beam could cause 
discrepancies in bone measurements that could be easily 
observed between serial radiographs. 10 Consequently, some 
researchers could have reported measurement discrepan- 
cies as bone loss or gain. 11 However, our preliminary study 
provides a foundation for future work and gives technical 
guidance for a prospective study on the scope of assess- 
ment of the trabecular bone in the periapical region of 
dental implants. A majority of similar studies have focus- 
ed on the marginal bone height and the osseointegration 
phase. 12 " 14 

Our results showed a significant difference between 
some parameters in the baseline and follow-up sets of 
radiographs in terms of the increase in the value of the 
quantifiable changes in the trabecular bone structure around 
the Frialit-2 implant after being put into clinical function. 
This can be masked by few data throughout the individual 
cases, meaning that although some cases showed an incre- 
ase in the follow-up values while others did not, there was 
an increase in the overall results. However, it can be argued 
that this difference could be related to the lack of perfect 
standardization measures in the procedure for taking radio- 
graphs. 

When the data from all of the patients were considered, 
there was no significant difference between the baseline 
and the follow-up radiographs for almost all the bone para- 
meters except the bone area (B Ar) and the mean distance 
between nodes (mean Nd-Nd). Thus, the null hypothesis 
was rejected for these two parameters and accepted for 
the others, indicating that there were no significant quanti- 
fiable changes in the trabecular bone structure around the 
endosseous implants after being put into clinical function 
for almost all of the parameters. Regarding the marrow 
spaces, the null hypothesis was accepted for all the para- 
meters except for bone area (BAr) and trabecular length 
(TbLe). What do these results mean? Seemingly, the most 
convincing result was the increase in bone area in the ROIs 
and the concordant fall in the marrow area. Thus, it appear- 
ed that the bone around the implant underwent remodel- 
ing that resulted in a net increase in the bone after implant 
placement. Therefore, the effect of successful implants on 
the bone did not appear to be limited to the immediate 
implant/bone contact area. If this result was correct, then 
the likely explanation was the effect of stresses and strains 
of implant function on the bone. However, this could be 
related to the scientific fact that the bone is a complex 
organ with cells specifically devoted to a continuous pro- 



cess of remodeling. According to Wolff's law, bone has 
the ability to alter its strength through the production of 
collagen fibers and increased deposition of minerals as a 
protective measure. Therefore, disuse of bone leads to atro- 
phy, whereas increased stresses lead to bone deposition. 
Another probable effect of mechanical stress is the incre- 
ased production of the calcitonin hormone that inhibits 
bone resorption. This information could be of great signi- 
ficance in explaining the increase in trabecular bone area 
(BAr) around Frialit-2 implants and relating the success 
of the implant to the equal distribution of stress around it. 

Our cases involved implant placements into the edentu- 
lous areas of varying durations. Therefore, it was likely 
that a degree of bone loss in the implant site space occurred 
during the intervening period. Thus, we postulated that 
implant placement and subsequent functional stresses 
caused the bone deposition detected by our analysis. It 
would be interesting to see whether this observation would 
be identifiable in cases of immediate implant therapy, 
where there would be no intervening time allowing for 
bone atrophy. 

In conclusion, several methods should be used in paral- 
lel to measure the trabecular architecture, as a single tech- 
nique may focus on a given aspect and neglect the others. 
In the context of our study, identification of alternative 
techniques was more difficult. The small scale and fine 
details on dental intraoral images implied that other sug- 
gested alternative techniques (such as peripheral microto- 
mographs) are probably not applicable. Rather, it would 
be possible to perform a subjective image analysis of the 
architectural or densitometric measurement. 
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